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ABSTRACT 

Context. Solar magnetic fields are regularly extrapolated into the corona starting from photospheric magnetic measurements that can 
suffer from significant uncertainties. 

Aims. Here we study how inaccuracies introduced into the maps of the photospheric magnetic vector from the inversion of ideal and 
noisy Stokes parameters influence the extrapolation of nonlinear force-free magnetic fields. 

Methods. We compute nonlinear force-free magnetic fields based on simulated vector magnetograms, which have been produced by 
the inversion of Stokes profiles, computed from a 3-D radiation MHD simulation snapshot. These extrapolations are compared with 
extrapolations starting directly from the field in the MHD simulations, which is our reference. We investigate how line formation 
and instrumental effects such as noise, limited spatial resolution and the effect of employing a filter instrument influence the resulting 
magnetic field structure. The comparison is done qualitatively by visual inspection of the magnetic field distribution and quantitatively 
by different metrics. 

Results. The reconstructed field is most accurate if ideal Stokes data are inverted and becomes less accurate if instrumental effects and 
noise are included. The results demonstrate that the non-linear force-free field extrapolation method tested here is relatively insensitive 
to the effects of noise in measured polarization spectra at levels consistent with present-day instruments. 

Conclusions. Our results show that we can reconstruct the coronal magnetic field as a nonlinear force-free field from realistic 
photospheric measurements with an accuracy of a few percent, at least in the absence of sunspots. 

Key words. Sun: corona - Sun: magnetic fields - Sun: photosphere 



1. Introduction 

Except for a few individual cases , e.g., in newly developed ac- 
tive regions (ISolankietalJl2003h . we cannot measure the full 
magnetic vector in the solar corona at high resolution directly 
and we ha ve to rely on extrapolation s from photospheric mea- 
surements. |\ViegejmamretiO^^ direct mag- 
netic field measurements by ISolanki et al.l d2003l) with extrap- 
olations from the photosphere. This work revealed the impor- 
tance of field aligned electric currents for an accurate mag- 
netic field reconstruction and the need for photospheric vector 
magnetograms as boundary data. The photospheric data suf- 
fer from a number of inadequacies, however, whose influence 
on the quality of the extrapolations needs to be studied in de- 
tail. Firstly, the magnetic fi eld in the photosphe re is far from 
being force-free, (see, e.g., iMetcalf et al.l Il995l) . which leads 
to significant errors if used directly for a force-free magnetic 
field extrapolation (see IMetcalf et al.l I2008L for details). These 
can, however, be g reatly r educed by appropriate preprocessing 
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dWiegelmann et al]|2006bl) . Another well known problem is that 
the noise level of the magnetic field transverse to the line of sight 
is typically more than one order of magnitude higher than for 
the longitudinal component deduced from the Stokes parame- 
ters of Zeeman split spectral lines. An additional complication 
arises for vector magnetic field measurements made far away 
from the disk center, where the vertical and line-of-sight field 
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are far apart (see lGarv & Hag vardl 19901 for details). In this case 
one has to use also the transvers e field to derive the vertical field. 
IVenkatakrishnan & Gary! d 1989b showed that the increased noise 
by using the transverse field is tolerable for heliocentric dis- 
tances of less than 50°. For our investigations we assume that the 
observations are made close to the disk center, where the vertical 
magnetic field is almost identical with the line-of-sight field and 
the horizontal component is almost identical with the transverse 
field. Consequently our results are not directly applicable to re- 
gions observed for heliocentric distances of larger than about 
50°, where the noise in the transverse field influences signifi- 
cantly the accuracy of the vertical field. Furthermore, high reso- 
lution vector magnetographs such as Hinode/SOT have a limited 
field of view, which often does not allow deriving the horizon- 
tal magnetic field in an entire active region. The influence of 
these effects on the accuracy of non-line ar force-free extrapo - 
lations have recently been investigated in iDeRosa et al.l d2009l) . 
Less often considered are effects introduced by the fact that the 
extraction of the field from the Stokes parameters is intrinsi- 
cally uncertain. E.g., the measured Stokes profiles are formed 
in a highly dynamic atmosphere with a complex thermal and 
magnetic structure, while the normally applied inversion meth- 
ods impose r ather restrictive a ssumptions (e.g. Milne-Eddington 
atmosphere. lAuer et al.l 1 1 977b . In addition, due to the spatially 
fluctuating height of formation of the lines, the obtained mag- 
netic vector refers to different heights in the atmosphere at dif- 
ferent locations. For extrapolation methods one has to assume a 
single, often planar, height as the boundary condition. Finally, 
instrumental limitations impose restrictions, such as limited spa- 
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tial and spectral resolution, spectral sampling (in particular for 
filter instruments) and noise. 

Here we investigate how the extrapolated coronal magnetic 
field is influenced by noise and other instrumental artifacts (spa- 
tial resolution, limited spectral sampling), as well as the general 
limitations of the inversion of the measured Stokes profiles. To 
do so we use the results of 3D radiation MHD simulations. We 
compute synthetic lines from the data cubes of the relevant phys- 
ical quantities, add noise and apply the influence of typical in- 
strument parameters. Finally we invert these artificial measure- 
ments to derive synthetic vector magnetograms, which are then 
used as boundary conditions for a nonlinear force-free magnetic 
field extrapolation. We compare the reconstruction from ideal 
data, taken directly from the MHD simulations, with extrapola- 
tions starting from data containing instrumental effects and noise 
of different levels of severity. Our aim is to investigate how dif- 
ferent instrumental effects and noise influence a nonlinear force- 
free magnetic field extrapolation. Besides testing generally the 
influence of reduced spatial resolution of the photospheric mag- 
netic field data, we also consider to what extend the properties 
of specific high resolution space instruments affect the quality 
of the extrapolations. The two instruments we consider are the 
Spectro-Polarimeter (SP) o n the Solar Opt i cal Telescope (SOT) 
on the Hinode spacecraft dShimizul 12004 iTsuneta et all 120081) 
and the Polarimetric and Helioseismic Imager (PHI) to fly on 
the Solar Orbiter mission. Instruments with lower spatial resolu- 
tion, such as the Helioseismic Magnetic Imager (HMI) were not 
considered, since the number of pixels across the MHD simula- 
tion box are then rather small, making a meaningful test more 
improbable. 



2. Method 

2. 1 . Setup of the test case 

We start with 3-D r adiation MHD simu lations resulting from the 
MURAM code (see lVogler et al.ll2005l for details). The particu- 
lar snapshot considered here has been taken from a bipolar run 
and harbours equal amounts of magnetic flux of both polari- 
ties. The configuration has an average field strength of 150 G 
at the spatially averaged continuum optical depth (tjooo) = 1 
at 5000 A. The horizontal dimensions of the simulation domain 
are 6 x 6Mm and 1.4Mm in the vertical direction. The origi- 
nal resolution is 20km in the horizontal directions and 14km in 
the vertical direc t ion. Sim ilar snapshots have also been used by 
iKhomenko et al.l (l2005bl) for the investigation of magnetocon- 
vection of mixed-polarity q uiet-Sun regions and com pared with 
measured Stokes profiles by IKhomenko et all (12005 al) . 

Due to the limited height-range of the MHD simulation we 
first extrapolate the field from the magnetic vector obtained 
from the simulations at a fixed geometric height (see Fig. Q]). 
This height, Z re f, must be chosen with care since the extrapola- 
tion starting from the magnetic field vector at Z re f is employed 
as a reference against which all others are compared. In or- 
der to avoid introducing a bias it must correspond to a suitably 
weighted average height of formation of the spectral line used 
for inversion, e.g. Fe I 6173.3 A. Since different regions in the 
MHD-simulation may have different formation height ranges, 
we have chosen a criterion to select the reference height which is 
obtained by finding the highest correlation between the magnetic 
field strength taken directly from the MHD simulation and from 
inversion of the Stokes profiles obtained under ideal conditions 
(very high spectral and spatial resolutions and no noise). 



This extrapolation is used as reference, with which other ex- 
trapolations are compared. These later ones use magnetic vector 
maps that are obtained from the inversion of synthetic Stokes 
profiles. 

Two spectral lines are considered, the very widely used Fe I 
6302.5 A (Lande factory = 2.5) line, employed by the Advanced 
Solar Polarimeter and the spectropolarimeter on the Hinode 
spacecraft and the Fe I 6173.3 A (g = 2.5), which has been 
selected for the Helioseismic and Magnetic Imager (HMI) on 
the Solar Dynamics Observatory (SDO) and for the Polarimetric 
and Helioseismic Imager on the Solar Orbiter mission (SO/PHI). 
For both Z re f was found to lie approximately at 150 km above 

<T"5000) = 1- 

Next, using the STOPRO code dSolanldll987tlFrutiger etail 

|2000|) . we compute the Stokes parameters of these widely used 
spectral lines. We restrict ourselves to considering the situation 
at solar disk center. These synthetic data are then either directly 
inverted, or after degrading the data by introducing noise, lower- 
ing the spatial and/or the spectral resolution, etc. A list of all the 
considered cases is presented in Table 1 and discussed in Sect. 
3. The inversion is carried out for each spectral line individually. 
The inverted magnetic field vector is then taken as the starting 
point for a non-linear force-free extrapolation. 

For comparison we also compute a potential field from the 
results of the MHD-simulations at Z re f , which requires only the 
vertical component of the photospheric magnetic field. In ad- 
dition, we compute nonlinear force-free extrapolations with the 
horizontal spatial resolution reduced by a factor of two and four, 
respectively, compared to the original MHD model. 



2.2. Description of the inversion code 

The inversion code HeLIx developed by lLagg et all d2004l) is 
based on fitting the observed Stokes profiles with synthetic ones 
obtained from an analytic solution of t he Unno-Rachkovsky 
equations (Unno 1956; Rach kowskvll967|) in a Mil ne-Eddington 
atmosphere (see e.g. lLandi deglTnnocentilll992l) . These syn- 
thetic profiles are functions of the magnetic field strength \B\, 
its inclination and azimuth (with respect to the line of sight), the 
vertical velocity, the Doppler width, the damping constant, the 
ratio of the line center to the continuum opacity, and the slope of 
the source function with respect to optical depth. We assume the 
filling factor to be unity. Any possible magnetic fine structure 
within one pixel consisting of a field-free and a magnetic area is 
smeared out, and the inversion returns the magnetic field vector 
averaged over the considered pixel. Since the extrapolations are 
performed on the same pixel scale as the inversions, the extrap- 
olations must be fed with the pixel-averaged magnetic field vec- 
tor. Setting the filling factor to unity already during the inversion 
process enhances its robustness and reliability. The atmospheric 
parameters that ensure the minimum of the meri t function are 
obtain ed using a very reliable genetic algorithm dCharbonneaul 
119951) . The genetic algorithm has been extensively tested with 
synthetic spectra from MHD simulations and the result com- 
pared with response-function-averaged physical parameters (e.g. 
magnetic field strength, inclination, azimuth, line-of-sight veloc- 
ity). The results of the test indicate that the genetic algorithm 
retrieves the global minimum of the merit function with high re- 
liability. The parameter sets are chosen to, e.g., mimic the effects 
of the Hinode/SOT and the future SO/PHI instruments, respec- 
tively. 
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Fig. 1. B x ,By,B- maps at a fixed reference height Z le f approximately 150 km above the average < T5000 >= 1- The top panel 
contains the original maps from the MHD data, which are used as a reference. The second panel shows the results of an inversion 
with Filter and noise. The third and bottom panel correspond to simulated inversions from Hinode/SP and SO/PHI, respectively. For 
a quantitative comparison see Table [TJ 



2.3. Hinode SOT Spectro-polarimeter 

The spectropolarimeter (SP; iLites et aTl 120071) is part of the 
focal-plane package of the 50cm Solar Optical Telescope (SOT) 
onboard the Hinode spacecraft. It observes the line pair Fe I 
6302.5 and 6301.5 A. Here we restrict ourselves to 6302.5 A. 
The spatial resolution at the diffraction limit of the telescope's 
primary mirror is about 0.32" at 6302.5 A, which corresponds to 
230 km on the Sun. The size of a detector pixel corresponds to 
approximately 110 km on the Sun in the spatial direction. The 
spectral resolution is 30 mA and the spectral sampling is 21 mA. 
We used Gaussians for spectral and spatial smearing. The noise 
has been added to the Stokes profiles as photon noise Fxl/ y/F-, 



where F is a white noise with a Gaussian distribution and I c is 
the continuum intensity. The chosen standard deviation of F was 
10 3 , which corresponds to the typical noise level for modern 
spectropolarimetric observations (e.g. Hinode/SP data). Noise is 
included in a similar way wherever we mention that noise has 
been added to the data. More subtle instr umental effects such a s 
scattered light, or a slight defocus (e.g. iDanilovic et aT1l2008l) . 
are not considered. 



2.4. SO/PHI magnetograph 

The Solar Orbiter Polarimetric and Helioseismic Imager 
(SO/PHI), a vector magnetograph, will be one of the main instru- 
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merits on the ESA-NASA Solar Orbiter mission. Of its two tele- 
scopes, the High Resolution Telescope (HRT) is of primary im- 
portance for the present investigation (due to the small horizon- 
tal extent of the MHD simulations). The spectral line chosen for 
SO/PHI, Fe I 6173.3 A, combines high Zeeman sensitivity with 
spectral purity, needed for simultaneous vector magnetic field 
and helioseismology studies. We describe the point spread func- 
tion (PSF) by a Gaussian with FWHM= 150 km. The obtained 
arrays of Stokes parameters are rebinned to a spatial pixel size 
of 80km. We then convolve the Stokes profiles in the spectral di- 
mension with a Fabry-Perot type filter with FWHM = 100 mA. 
Since SO/PHI will be a filtergraphic instrument, we decrease the 
number of spectral samples per line by taking 5 positions in the 
line and one in the continuum at the positions (from line center 
at rest): -0.3A, -0.15A, -0.075 A, A, +0.075 A, +0.15 A. At 
that stage we add noise and perform the inversion of the Stokes 
profiles. 

2.5. Extrapolation of the vector magnetogram into the 
atmosphere 

To compute a 3D-nonlinear force-free magnetic field from the 
result of the HeLIx inversion code we carry out the following 
steps: 

- If needed, transform B, 9 and <p, which are output from the 
inversion code, to B x , B y , B z on the photosphere, which re- 
quires a resolution of the 180° ambiguity in <p. (See section 

ESU) 

- Preprocess the vector magnetogram (B x , B y , B z ), assuming 
that it refers to the same geometric height at every spatial 
pixel. (See Sect. l2.6.TI ) 

- Compute a nonlinear force-free coronal magnetic field from 
the preprocessed vector magnetogram. (See Sect. 12.6.21 ) 

- Compare the result with the reference field. 

We explain these steps in the following. 
2.5.1. Removal of the 180° ambiguity. 

For the purpose of the present investigation we have chosen to 
remove the 180° ambiguity by minimizing the angle to the exact 
solution. This possibility does, of course, not exist for real data 
and other methods for the ambiguity inversion would need to 
be tried. The performance of different ambiguity removal tech- 
niques has been studied with synthetic data by iMetcalf et alj 
d2006l) . They found that the best available technique managed 
to get 100% of the points correctly. Recently the influence of 
noise and spatial resolution on the quality of the di fferent am- 
biguity removal techniques has been investigated by iLeka et aT] 
(2009). We have therefore not considered specific ambiguity re- 
moval techniques. An investigation of their efficiency and influ- 
ence is outside the scope of this paper. 

2.6. Effects in the photosphere 

In Tableland in figure[T}ve investigate how the different instru- 
ment effects and noise influence the vector field in the photo- 
sphere. The first line corresponds to the MHD reference field. 
The field has an average electric current density of (\J Z \) = 
0.059A/m 2 and a vertical magnetic field strength of (\B Z \) = 
103.5G. Positive and negative values of these quantities are bal- 
anced. In the table we compute the correlation relative to this 
reference case for the horizontal fields Cb x , Cb v , the vertical 



Table 1. Influence on the photospheric field for different cases. 



Case studied C flA C By C Bz C vec < \J-\ > < \B Z \ > 



Hinode/SP 














no noise 


0.80 


0.77 


0.84 


0.82 


0.047 


79.4 


1(T 3 noise 


0.76 


0.71 


0.84 


0.81 


0.067 


80.7 


SO/PHI 














no noise 


0.79 


0.68 


0.85 


0.82 


0.078 


84.5 


1CT 3 noise 


0.67 


0.53 


0.85 


0.73 


0.184 


85.2 



field Cb z and the 2D-vector correlation in the photosphere C vec 
and provide the average absolute values {\B Z \) and {\J Z \). For all 
full spatial resolution cases (upper part of the table) the correla- 
tion for the different cases is 0.97 for the vertical magnetic field 
strength and the average field strength is underestimated by a 
few percent. Noise and instrument effects seem to have a relative 
small effect on the vertical field. The effect of a reduced resolu- 
tion (lower part of the figure, Hinode/SP and SO/PHI cases) on 
the horizontal photospheric magnetic fields B x , B y and the de- 
rived vertical current density J z is significantly higher. While for 
full profiles the correlation in the horizontal fields is in the range 
of 0.85 - 0.88, the combined effect of filter and noise reduces 
the correlation to only 0.70 and spurious, non-physical electric 
currents, which results in an overestimation of J z by a factor of 
about 3. 

2.6.1. Preprocessing 

The magnetic field in the photosphere is not necessarily force- 
free ( because of the finite B plasma in the photosphere, see lGarvl 
|2001|) and the horizontal components (B x an d B v ) of cur rent vec- 
tor magnetographs have large uncertainties. lAlvl dl989l) defined 
a number of integral relations to evaluate if a measured photo- 
spheric vector magnetogram is consistent with the assumption 
of a force-free field. These integral relations (numerator in Eq. 
Q3 have been used to define a dimensionless parameter ef orce as 

| JT B X B Z dxdy\ + | J s B y B z dxdy\ + | J S (B X + B y ) - B\ dxdy\ 

J S (B 2 X + Bj + B^)dxdy (1) 

and force-free extrapolation codes require ff orce « 1 on the 
boundary. (One gets £f orce = if the integral relations are ful- 
filled exactly). 

For the synthetic magnetograms investig ated in table [T| we 
find gf o rce = 0.57 ± 0.19 in the photosphere. IWiegelmann et al.l 
d2006b1) developed a preprocessing procedure to drive the ob- 
served non force-free data towards boundary conditions suit- 
able for a force-free extrapolation. As a result of the prepro- 
cessing we get a boundary-data set which is consistent with 
the assumption of a force-free magnetic field. After applying 
the preprocessing-routine we get force-free consistent boundary 
conditions with ff orce = (2.6 + 0.2) • 10~ 4 . The preprocessing af- 
fects mainly the horizontal magnetic field components. The cor- 



Reference 1 1 1 

Inversion of synthetic profiles (tests) 
Full resolution, no noise and full profiles 

6302.5 A 0.89 0.88 0.97 

6173.3 A 0.89 0.88 0.97 

Full resolution, with noise and full profiles 

6302.5 A 0.87 0.86 0.97 

6173.3 A 0.88 0.88 0.97 

Full resolution, with Filters (5 + I A) values 
6173.3 A no noise 0.86 0.85 0.97 
6173.3 A w. noise 0.72 0.70 0.97 0.90 0.213 



1 


0.059 


103.5 


0.96 


0.065 


96.7 


0.96 


0.065 


95.9 


0.95 


0.086 


97.6 


0.95 


0.070 


96.5 


0.95 


0.084 


96.0 


0.90 


0.213 


96.7 



Wiegelmann et al.: Force-free modelling: influence of spectral lines 



5 



Table 2. Various figures of merit of the 3D reconstruction for 
different cases. 



Vector correlation: 



Case studied 


r 

^vec 


e 


ClOOkm 


z ■ 

<- 400km 


L 800km 


Reference 


1 


1 


1 


t 
1 


t 
1 


Potential 


0.68 


0.62 


0.72 


U.o4 


n Q7 
u.y / 


MHD cases with reduced resolution 






Pixel Size 40km 


0.99 


0.93 


0.99 


1.00 


1.00 


Pixel Size 80km 


0.93 


0.79 


0.94 


0.96 


0.99 


Inversion of synthetic profiles (tests) 






Full resolution, no noise and full profiles 






6302.5 A 


0.92 


1.06 


0.93 


0.98 


1.00 


6173.3 A 


0.93 


0.99 


0.93 


0.98 


0.99 


Full resolution, with noise and full profiles 






6302.5 A 


0.93 


1.04 


0.93 


0.98 


1.00 


6173.3 A 


0.93 


1.05 


0.93 


U.Vo 


Vj.yy 


Full resolution, with Filters (5 + 1/t) values 






6173.3 A no noise 


0.92 


1.06 


0.92 


0.98 


0.99 


6173.3 A w. noise 


0.91 


1.10 


0.89 


0.97 


0.99 


Hinode/SP 












no noise 


0.85 


0.66 


0.85 


0.96 


0.99 


10~ 3 noise 


0.84 


0.62 


0.84 


0.95 


0.99 


SO/PHI 












no noise 


0.84 


0.75 


0.85 


0.96 


0.99 


10~ 3 noise 


0.81 


0.76 


0.78 


0.90 


0.98 



relation of original and preprocessed field for the investigated 
cases is Cg z = 0.99+0.006 for the vertical component B z . For the 
horizontal components we find a correlation of Cb.k = 0.91 ±0.05 
and C By = 0.90 + 0.05, respectively. 

2.6.2. Extrapolation of nonlinear force-free fields. 

Force-free coronal magnetic fields have to obey the equations 

(VxB)xB = 0, (2) 
V B = 0. (3) 

We define the functional 



L 



= f w [B" 2 |(VxB)xB| 2 + |V-B| 2 ] d 3 x, 



(4) 



where w is a weighting function. It is obvious that (for w > 0) the 
force-free equations (|2] |3]l are fulfilled when L is equal to zero 
dWheatland et a l. 2000). We min imize the functional © numer- 
ically as explained in detail by IWiegelmannl d2004l) . The pro- 
gram is written in C and has be en par allelized with O p enMP . 
IWiegelmann & Neukirchl (120031) and ISchrijver et all (120061) 
tested the program with ex act nonlinear force - free equ ilibria of 
lLow & Loul dl990l) . while IWiegelmann et all d2006al) tested it 
with a nother exact equilibrium developed by iTitov & Demoulinl 
d 19991) . The code has been applied to vector magnetograph 
data from the German V acuum Tower Telescope (VTT) by 
IWiegelmann et al] (2005a b) and to data from the Solar Flare 
Telescope (SFT) by IWiegelmann et al.l (l2006bl) . Here we use 
an updated version of the optimization approach, includ- 
ing a m ulti-scale approach which has been described and 
tested by iMetcalf et at] d2008l) and applied to Hinode data by 
dSchrijver et alj|2008l) . 

2. 7. Figures of merit 

ISchrijver et al.l d2006l) introduced several figures of merit to 
compare the results of magnetic field extrapolation codes (a 3D- 
vector field b) with a reference solution B. 



C vec = ^Bi-bi/ ^|Bi| 2 J]lb 



1/2 



(5) 



where i corresponds to all grid points in the entire 3D com- 
putational box. 

- Total magnetic energy of the reconstructed field b normal- 
ized by the energy of the reference field B: 



E/IBil 



(6) 



- We also compute the linear Pearson correlation C of the total 
magnetic field strength |B| at the heights 100km, 400km and 
800km above the reference height, Z re f, respectively. 

Two vector fields agree perfectly if C vec , e and the Pearson cor- 
relation coefficients are unity. 

3. Results 

Table [2] contains the different quantitative measures which com- 
pare the different reconstructed field from synthetic observations 
with the reference field (extrapolations from ideal data). Column 
2 contains the vector correlation C vec , column 3 the normalized 
magnetic energy e, and columns 4-6 the linear Pearson corre- 
lation C of the total magnetic field strength \B\ at the heights 
100km, 400km and 800km above the reference height, Z re /, re- 
spectively. 

Figs|2]and[3]show contourplots of the magnetic field strength 
|B| at 100 km and 400 km above Z re f respectively. 

3.1. MHD cases 

In the first 4 images of Figs [2] and [3] we compare the nonlinear 
force-free 3D magnetic field reconstructed from ideal data (ex- 
tracted from the MHD simulations at Z re f, called Reference in 
Table [2] and in Figs. [2] [3]) with a potential field extrapolation 
also starting from ideal data (Potential) and nonlinear force- 
free computations starting from magnetic field maps obtained 
from the inversion of synthetic Stokes profiles and MHD cases. 
Potential fields are often calculated in addition to NLFFF, be- 
cause they contain the minimum energy for given vertical bound- 
ary conditions and the free energy of a NLFF-field above that of 
a potential field has relevance for coronal eruptions. 

As seen in Table [2] and in Figs. [2] and [3] the computations 
with reduced spatial resolution of 40 km pixels show almost 
perfect agreement, while 80 km pixels still provide an excel- 
lent correspondence. The correlation with the reference field is 
better higher in the atmosphere, which is not surprising since 
small-scale field structures at the lower boundary do not propa- 
gate very high and the larger scales are less affected by the bin- 
ning to larger horizontal pixels. The magnetic energy is under- 
estimated because small scale fields and currents at low heights, 
which contribute significantly to the total magnetic energy, are 
not resolved here. 



3.2. Dependence on spectral line, noise and spectral 
sampling 

The remaining images in Figs. [2] and [3] and the remaining en- 
tries in Table [2] apply to extrapolations starting from magnetic 
maps obtained from inversions of synthetic spectral lines. The 
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Reference 



Pixel 40km 



Pixel 80km 





Potential 



6173, no noise 



6173, Filter, no noise 




6173, Filter, noise 



Hinode/SP 



SO/PHI 





Bl/Gauss 



Height 100 km 



Fig. 2. |B| for the reference field and extrapolated from inversion of synthetic profiles. The plots display |B| at a height of 100km 
above the reference height (See Sect. 12. II for a definition). The figure shows, from top-left to bottom-right, the ideal reference 
solution (indicated as Reference); extrapolations starting from the MHD output after reducing its spatial resolution (Pixel 40km 
and Pixel 80km, respectively); a potential field reconstruction (Potential); a reconstruction from a full resolution spectropolarimeter 
in the 6173.3 A line without noise and inverting the full profiles (6173, no noise); reconstructions starting from inversions of full 
resolution filter magnetograms without and with noise (Filter, no noise and Filter, noise, respectively); and finally two cases with 
specifications adapted to those of two space-borne instruments, the spectropolarimeter on Hinode (Hinode/SP) and the Polarimetric 
and Helioseismic Imager on Solar Orbiter (SO/PHI), x and y axis are in Mm. 



first 3 of these images (marked 6173 no noise, Filter no noise, 
Filter noise in Figs.[2]and[3]i and the whole middle part of Table 
|2]are test cases for which we move step-by-step away from ideal 
conditions (i.e. B obtained directly from the MHD simulations) 
towards more realism, i.e. introducing various instrumental ef- 
fects. 



The first such step is to invert the synthetic line profiles di- 
rectly, without any further manipulation (ideal spectropolarime- 
ter). This step introduces uncertainties due to the inversion pro- 
cess per se (e.g. by the fact that asymmetric Stokes profiles are 
fit assuming purely symmetric or antisymmetric ones), the fact 
that the line samples a range of heights, and the dependence of 
the line formation height on the type of solar feature. 
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Fig. 3. Same as Fig. [2] but AQQkm above the reference height. 



The inversion of Stokes spectra expected from an ideal spec- 
tropolarimeter leads to an agreement with the reference within 
a few percent for the vector correlation and we can estimate the 
magnetic energy with an accuracy of one percent for the 6173.3 
A line. The magnetic energy is an important quantity because it 
tells us how much free energy is maximally available for erup- 
tive phenomena like flares and coronal mass ejections. The high 
accuracy achieved in this case is particularly encouraging since it 
shows that applying a simple Milne-Eddington inversion to the 
often highly asymmetric profiles (sometimes showing multiple 
lobes in Stokes V) gives sufficiently accurate results and that the 
fluctuating height at which the magnetic field is obtained does 
not significantly influence the extrapolations (note that the situ- 
ation may be different in a sunspot with its rather deep Wilson 



Height 400 km 



depression). This result implies that most of the inaccuracies in 
the extrapolations are due to limitations in the instrumentation 
and not because it is in principle not possible to extract the in- 
formation from the observations. 



As Table [2] shows there is little to choose between the two 
g = 2.5 Zeeman triplets /16302.5 A line and ,16173.3 A. The 
field extrapolated from the magnetic field maps derived from ei- 
ther of them correlate very well with the reference field. Using 
the /16173.3 A lines leads to a slightly better estimate of the mag- 
netic energy, but all in all the necessity of using the Zeemann 
effect in a spectral line by itself only leads to errors of a few per- 
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centQ. The addition of noise at a level of 10~ 3 of the continuum 
intensity I c , which is typical of modern spectropolarimetric ob- 
servations, has only a small effect. Note that by adding an equal 
amount of photon noise to all Stokes parameters we are produc- 
ing a much lower S/N ratio in the linearly polarized Q and U 
profiles than in Stokes V due to the smaller signal in the former. 
The influence of noise of a given amplitude also depends heav- 
ily on the magnetic flux in the box for which we extrapolate. 
With (B) 150 G the chosen snapshot corresponds to an aver- 
age plage region. We expect that the same amount of noise will 
have a considerably stronger effect when applied to the weaker 
Stokes profiles present in the quiet sun. 

It is of particular interest if measurements in limited wave- 
length bands and at reduced wavelength resolution, typical of 
filter polarimeters (filter magnetographs), are acceptable for ex- 
trapolations. One advantage of such instruments is that they al- 
low time series of a whole region to be recorded at high cadence. 
Of particular relevance for magnetic extrapolations is the fact 
that filter polarimeters record the Stokes vector over the full field 
of view. This overcomes the main shortcoming of spectropo- 
larimeters, namely that they need to scan a region step by step, 
so that by the time the second footpoint of a loop is scanned the 
first may have evolved considerably. These advantages of filter 
instruments come at the price of a reduced spectral resolution 
and limited spectral sampling. An extensive series of tests by 
one of us (L. Yelles) using various MHD simulation snapshots 
has shown that observations at 5 wavelength points in the line 
plus one at the continuum should be adequate to obtain the mag- 
netic field vector reliably. 

The computations carried out here suggest that this is also 
true for the magnetic field extrapolated from vector magne- 
tograms obtained from filter instruments. According to Table [2] 
and Figs.|2]and|3]the application of 100 mA broad filters to 5 lo- 
cations in the Fe 1 6173.3 A line and additionally to a nearby con- 
tinuum position gives an extrapolated nonlinear force-free field 
that differs only slightly from the results obtained with the full 
line profile. The magnetic energy is overestimated by at most 6% 
for a filter width < 100 mA and at least 5 sampling points in the 
line and one in the nearby continuum. 

Noise has a somewhat larger effect on filter polarimeter mea- 
surements than on spectropolarimetric ones, as can be judged 
from Table [2] In particular, the magnetic energy is affected, be- 
ing more than 2x less accurate than for spectropolarimetric mea- 
surements. If we take the effect of photometric noise of 10 3 / c 
in all Stokes parameters into account we overestimate the mag- 
netic energy by 10%. The somewhat higher magnetic energy for 
these cases is probably a result of stronger currents in the pho- 
tosphere, which are created by a less accurate computation of 
the horizontal photospheric field during the inversion. The total 
vertical magnetic field J \B Z \ dxdy in the photosphere is underes- 
timated by 7% after inversion of a set of filter images (irrespec- 
tive of whether noise is applied or not). The total vertical current 

f\ltx ~ IJv'l dxdy is overestimated by 42% and 262% without 
and with noise, respectively. Most of these spurious currents pro- 
duced by spectral line inversions fluctuate on very small scales 
so that they either are noise, or behave like that. Consequently, 
they are not transported into the corona. The preprocessing rou- 
tine takes care of this problem and most of the spurious currents 



1 The Z rcf was set to the average formation height for the /16173.3 A 
line. This might explain the somewhat worse results for the /t6302.5 A 
line. 



vanish. After preprocessing the total current is overestimated by 
2% and 9% for inversions without and with noise, respectively. 

3.3. Hinode-like cases 

Finally, as discussed in this and the following subsection, we 
add more realism into the synthetic observations by employing 
parameters that are appropriate to high resolution instruments. 
In this section we consider the important case of the spectropo- 
larimeter on Hinode (see also Section |2~3l ). 

Taking a finite spatial resolution (pixel size 110 km) and 
spectral smearing into account naturally leads to less accurate re- 
sults (rows marked Hinode/SP in Table|2]i. We find that we can- 
not reconstruct the magnetic energy accurately for these cases, 
because small scale magnetic fields are not adequately resolved. 
We get, however, a reasonable estimate of the correct magnetic 
field in higher layers of the atmosphere. For the most involved 
case (Hinode-like spectropolarimeter with a noise of 10 3 / c and 
a pixel size of 1 10 km) we find an error of 16%, 5% and 1% at 
the heights 100, 400 and 800 km above Z re f, respectively. The 
limited resolution avoids accurate reconstructions of low-lying 
small scale features, but gets the higher lying field approximately 
correctly. At heights above 400 km even the most involved and 
noisy nonlinear force-free reconstruction considered here has an 
accuracy, which is three times better than with a potential field 
reconstruction starting from a perfectly known lower boundary. 

3.4. SO/PHI-like cases 

We consider instrumental effects appropriate to the PHI instru- 
ment on Solar Orbiter such as the finite spectral resolution and 
sampling by a Fabry-Perot interferometer, finite spatial resolu- 
tion of approximately 160 km (pixel size 80 km) on the Sun 
and photon noise at a level of 10~ 3 x -^=. Details are given in 

Sect. 12.41 The extrapolated magnetic field displays a very simi- 
lar spatial distribution as that based on the inversion of ideal line 
profiles. On the whole, the influence of degrading the spatial res- 
olution to a pixel size of 80 km (denoted as Filter +Noise in Figs. 
[2] [3]) introduces a similar level of inaccuracy in the extrapolated 
field as the uncertainties in B x , B y , B z at the lower boundary in- 
troduced by instrumental effects and the inversion of the line pro- 
files. But even for this most involved cases, with comparatively 
low resolution, instrument effects and noise, the agreement with 
the reference is better than a potential field computed from ideal 
data at the lower boundary. The better correspondence of e with 
the reference value than for Hinode is likely due to the somewhat 
higher spatial resolution expected for PHI. This suggests that for 
accurate estimates of the magnetic energy it is more important 
to achieve high spatial resolution than completely accurate line 
profiles. The field at 100 km above Z re f and partly at 400 km 
is better reproduced for Hinode/SP-like parameters than for the 
SO/PHI case. Obviously, for this it is better to have the full line 
profile (in particular when noise is included). 

4. Conclusions 

We have investigated how strongly inaccuracies in the lower 
magnetic boundary, introduced by measurement and analysis er- 
rors of the Stokes profiles influence the computation of nonlinear 
force-free coronal magnetic fields. 

We find that instrument effects and noise influence the hor- 
izontal component of the photospheric magnetic field vector 
stronger than the vertical field. We find that non-linear force- 
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Fig. 4. Scatter plot displaying how errors in the photosphere af- 
fect the accuracy of the 3D magnetic field reconstruction. Shown 
is the 3D vector correlation against the 2D vector correlation in 
the photosphere. The solid line corresponds to a linear fit and the 
dotted line corresponds to identical 2D and 3D correlations. 

free fields extrapolated from ideal data and from the inversion 
of Stokes profiles deviate more strong at low heights. In par- 
ticular, a limited spatial resolution influences the lowest lay- 
ers the strongest. Higher in the atmosphere we found a very 
good agreement (correlation better than 0.98) with extrapola- 
tions from ideal data. We find that for an accurate estimation of 
the magnetic energy a high spatial resolution is more important 
than a high spectral resolution. 

These basic findings apply to magnetic vector maps ob- 
tained from both, spectropolarimetric data, such as provided by 
Hinode/SP, and filter magnetographs, such as to be provided by 
SO/PHI or SDO/HMI. 

Finally, we would like to determine how errors in the pho- 
tosphere influence the quality of the 3D reconstruction. For this 
aim we show a scatter plot in Fig. [4] to compare the 2D-vector 
correlation in the photosphere (Table [TJ with the 3D vector cor- 
relation of the reconstructed coronal magnetic field (Table 
The solid line in Fig. [4] shows a linear fit to the data points and 
the dotted line corresponds to equal 2D and 3D correlations. As 
one can see from this figure the relation between the accuracy 
of the photospheric and coronal field is linear, but not identical. 
In particular, the correlation of the full 3D field drops off more 
slowly than the photospheric correlation, so that extrapolations 
based on a filter instrument with noise (which gives the least 
accurate photospheric field) are more accurate than the vector 
magnetograms they are based on. 

It would be interesting to repeat this study on larger spatial 
scales, reaching sizes typical of observed vector magnetograms, 
as soon as the corresponding radiative MHD-simulations be- 
come available (with a similar grid size as the simulations em- 
ployed here). There is also a need to investigate quiet Sun re- 
gions, coronal holes and active regions separately. In particular, 
our results may not apply when considering extrapolations start- 
ing from regions containing sunspots due to their much larger 
Wilson depression and rather different temperature structure, 
which influence the line formation height and the line forma- 



tion in general. We have also not considered subtle, but possibly 
important effects such as the evolution of the field during the 
scan of an active region by the slit of a spectropolarimeter, or 
the evolution of a line profile during the spectral scan of a filter 
instrument. 

One might also consider investigating the lower layers of the 
solar atmosphere, where the plasma is not force-free, in more 
detail and taking non-magnetic forces into consideration for the 
magnetic field extrap olation. A first step in this directi on has re- 
cently been taken by IWiegelmann & Neukirchl (I2006I) . who de- 
veloped a magnetohydrostatic extrapolation code. For an appli- 
cation to data we require, however, more information regarding 
the nature of the non-magnetic forces (e.g., pressure gradients 
and gravity). 
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